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Mid-latitude  Temperatures  at  87  km:  Results  From 
Multi-instrument  Fourier  Analysis 
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Abstract.  Using  a  novel  Fourier  fitting  method  we  com¬ 
bine  two  years  of  mid-latitude  temperature  measurements 
at  87  km  from  the  High  Resolution  Doppler  Imager,  the 
Colorado  State  University  lidar,  and  the  Peach  Mountain 
Interferometer.  After  accounting  for  calibration  bias,  sig¬ 
nificant  local-time  variations  on  the  order  of  10  K  were  ob¬ 
served.  Stationary  planetary  waves  with  amplitudes  up  to 
10  K  were  observed  during  winter,  with  weaker  wave  ampli¬ 
tudes  occurring  during  other  seasons.  Because  of  calibration 
biases  among  these  instruments,  we  could  estimate  the  an¬ 
nual  mean  temperature  to  no  better  than  193.5  =L  8.5  K. 

1.  Introduction 

State  parameters  in  the  mesosphere  and  lower  thermo¬ 
sphere  (henceforth  MLT)  are  difficult  to  estimate  because 
significant  observational  aliasing  can  occur  due  to  irregular 
and  incomplete  sampling  of  the  region  in  space  and  time 
[. Salby ,  1982;  Forbes  et  al,  1997].  These  problems  are  com¬ 
pounded  by  the  possibility  of  direct  observational  biases  spe¬ 
cific  to  a  given  measurement  technique.  While  advanced 
data  assimilation  methods  and  sufficient  data  exist  for  tro¬ 
pospheric  and  stratospheric  objective  analysis,  MLT  tem¬ 
perature  data  and  dynamical  models  are  currently  inade¬ 
quate  to  employ  such  techniques.  Two  notable  attempts, 
however,  to  combine  multi-instrument  satellite  and  ground- 
based  MLT  data  to  specify  the  spatiotemporal  variability 
of  the  region  are  the  MSISE-90  empirical  model  of  Hedin 
[1991]  and  recent  work  by  Leblanc  et  al.  [1999].  We 
present  a  methodology  for  making  statistical  estimates  of 
mid-latitude  temperature  fields  using  a  combination  of  satel¬ 
lite  and  ground-based  data.  Results  at  a  height  of  87  km 
are  presented  here  to  illustrate  the  method  and  its  potential 
future  use  with  larger  multi-instrument  data  sets. 

Figure  la  shows  the  combined  longitudinal  and  tempo¬ 
ral  sampling  pattern  of  the  temperature  field  between  41° 
=L  1°  N  at  87  km  over  a  two-year  period  from  the  High 
Resolution  Doppler  Interferometer  (HRDI)  [ Ortland  et  al , 
1998],  the  Colorado  State  University  (CSU)  lidar  [She  and 
von  Zahn ,  1998],  and  the  Peach  Mountain  Interferometer 
(PMI)  [Niciejewski  and  Killeen,  1995].  The  CSU  lidar  and 
PMI  are  located  at  40.6°  N,  105°  W  and  41°  N,  83°  W,  re¬ 
spectively.  The  OH(3-l)  rotational  temperatures  measured 
by  the  PMI  are  assumed  to  originate  from  a  layer  at  87  =L 
4  km  [She  and  Lowe ,  1998].  Additional  information  on  the 
mean  altitude  of  the  OH  nightglow  layer  and  its  spatiotem- 
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poral  variability  is  provided  by  Yee  et  al.  [1997]  and  Melo 
et  al.  [1999].  In  this  study  we  have  used  HRDI  level  3 AT 
version  11  data  which  have  a  vertical  resolution  of  3  km. 

Figure  lb  shows  the  corresponding  sampling  patterns  as 
a  function  of  day  number  and  local  time.  Gaps  are  due 
to  the  HRDI  yaw  cycle,  cloud  cover,  and  instrument  down¬ 
time.  The  HRDI  observations  are  daytime  measurements, 
while  the  ground-based  observations  are  performed  at  night. 
Figure  lb  illustrates  the  asynoptic  sampling  pattern  of  the 
HRDI  instrument,  which  precesses  through  all  local  times 
in  36  days.  Note  from  Figure  1  that  combining  data  from 
all  three  instruments  significantly  improves  the  spatial  and 
temporal  data  coverage  compared  to  consideration  of  data 
from  a  single  instrument. 

2.  Methodology 

Next  we  seek  to  combine  these  data  to  make  a  de¬ 
tailed  estimate  of  mid-latitude  temperature  variability  at 
87  km.  As  a  prototype  we  have  used  non-sequential,  error- 
weighted,  least-squares  estimation  without  dynamical  con¬ 
straints.  The  idea  is  to  perform  simultaneously  two  Fourier 
harmonic  surface  fits  of  the  space-time  distribution  of  mea¬ 
surements  (Figure  la)  and  the  time-time  distribution  (Fig¬ 
ure  lb).  The  inversion  method  combines  space-time  anal¬ 
ysis  to  capture  the  quasi-stationary  planetary  wave  field 
[e.g.  Salby ,  1982]  with  non-linear  time  series  analysis  to  cap¬ 
ture  the  migrating  tidal  field,  while  accounting  for  seasonal 
changes  in  the  mean.  An  important  zeroth-order  correction 
term  is  also  included  to  account  for  possible  instrument  bi¬ 
ases.  The  multi-dimensional  Fourier  model  is: 

N  S 

[ons  cos(«sA  +  ont)  +  bna  sin(«aA  +  <Tnt)] 

71  =  0  8=  —  S,8^0 

N  M 
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The  first  Fourier  expansion  is  a  space-time  model,  where 
t  is  time  in  days,  A  is  longitude  and  ana  and  bns  are  real 
Fourier  coefficients.  This  term  represents  the  inter-seasonal 
variations  in  the  longitudinal  structure  of  quasi-stationary 
planetary  waves.  The  spatial  and  temporal  frequencies  are 
Kg  =  s27t/360°  and  an  =  n27r/ATn  with  s  and  n  as  the 
respective  integer  harmonics.  The  windowing  period  A Tn 
and  total  orders  N ,  M,  and  S  can  be  chosen  to  suit  the  task 
at  hand.  In  the  second  Fourier  expansion,  r  is  local  time 
in  hours,  Cnm  and  dnm  are  real  Fourier  coefficients  for  inte¬ 
ger  harmonics  n  and  m,  respectively,  and  =  m2ir / ATm . 
This  term  represents  the  seasonal  amplitude  and  phase  mod¬ 
ulations  of  the  fundamental  westward  propagating  tidal  har¬ 
monics,  i.e.  the  diurnal  (1,1)  and  semi-diurnal  (2,2)  modes. 
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The  third  term  in  (l)  represents  the  seasonal  variation  of 
the  zonal  mean,  which  has  been  separated  from  the  first  two 
terms  to  avoid  redundant  harmonics.  The  last  term  is  in¬ 
troduced  to  account  for  zeroth  order  instrument  calibration 
biases,  where  a *.  is  an  indicator  function  which  takes  the 
value  of  unity  for  an  observation  from  the  kth  instrument 
(otherwise  it  is  zero).  It  is  used  to  tag  observations  from 
the  kth  instrument  for  multiplication  by  the  bias  correction 
factor  Pk  in  a  matrix  formulation  of  the  forward  model. 

With  a  total  of  J  measurements  from  K  instruments 
the  25  +  2 NS  +  2 M  +  2 NM  +  2 N  +  K  unknown  coeffi¬ 
cients  (dns i  bns;  Cnjn,  dnm^  fn^  Ph)  can  be  estimated  in  an 
err  or- weighted,  least-squares  sense  from  the  overdetermined 
system  of  J  linear  equations  using  Singular  Value  Decom¬ 
position  (SVD)  inverse  methods  [ Wunsch ,  1996].  Because 
all  harmonic  components  are  estimated  simultaneously  and 
are  resolved  by  the  data  (see  Figure  1),  aliasing  among  the 
mean,  quasi-st  at  ion  ary  planetary  waves,  and  tidal  compo¬ 
nents  is  effectively  eliminated. 

In  constructing  this  model  we  have  assumed  that  the  rel¬ 
ative  variations  within  a  given  data  subset  are  unbiased, 
but  that  time-independent  calibration  biases  are  possible 
in  each  data  set.  The  coefficients  pk  represent  the  aver¬ 
age  departure  of  observations  from  the  kth  instrument  from 
some  annual  mean.  To  facilitate  this  formulation  the  an¬ 
nual  mean  temperature  (estimated  here  as  the  sample  mean 
of  the  entire  data  set)  is  subtracted  from  the  data  before 
fitting.  Analysis  performed  using  synthetic  data  indicates 
that  adequate  information  exists  in  the  data  to  resolve  both 
the  diurnal  harmonic  and  the  bias  terms,  despite  the  func¬ 
tional  similarities.  If  time-dependent  instrument  biases  exist 
in  the  data  (for  example  from  seasonal  changes  in  the  height 
of  the  OH  nightglow  layer)  the  corresponding  effect  on  the 
state  estimates  will  be  weighted  inversely  by  the  statistical 
uncertainties  of  the  data  and  proportionately  by  the  rela¬ 
tive  influence  of  the  data  subset  on  each  model  parameter. 
In  future  analyses  of  upper- atmospheric  data  sets,  valuable 
insight  into  time-dependent  observational  biases  can  be  ob¬ 
tained  by  application  of  (1)  followed  by  a  statistical  analysis 
of  the  solution  residuals. 

A  number  of  choices  of  N,  5,  M,  ATn,  and  A Tm  are 
possible.  We  show  results  for  the  current  data  set  using 
N  =  5,  M  =  3, 5  =  5,  ATn  =  730  days,  and  A Tm  =  24 
hours.  In  the  current  inversion  the  solution  at  the  endpoints 
is  forced  to  be  periodic  so  that  the  estimates  for  January 
1,  1993  and  December  31,  1994  will  be  similar.  Over  the 
two-year  period  there  were  J  =  12685  observations;  7672 
observations  at  87  km  between  40°  and  42°  N  made  by 
HRDI,  2496  hourly- aver  aged  measurements  between  86.0 
and  88.0  km  made  by  the  CSU  lidar  (3  x  832),  and  2517 
hourly- averaged  measurements  made  by  the  PMI.  In  order 
to  account  properly  for  geophysical  variability  from  gravity 
waves  that  is  not  captured  by  (1)  and  more  evenly  weight 
the  three  data  sets,  the  reported  measurement  uncertainties 
were  adjusted  to  reflect  the  magnitude  of  geophysical  root 
mean  square  temperature  fluctuations.  While  maintaining 
the  relative  variations  within  each  data  set  the  reported  un¬ 
certainties  were  renormalized  to  a  mean  value  of  7  K. 

3.  Results 

The  sample  means  and  standard  deviations  of  the  indi¬ 
vidual  data  sets  over  the  two  year  period  are  186.6  =b  15.2  K 
(HRDI),  203.4  ±  18.8  K  (CSU)  and  205.8  ±  22.1  K  (PMI). 


Differences  in  these  means  can  be  interpreted  as  day- to- night 
temperature  changes,  instrument  bias  effects,  or  a  combina¬ 
tion  of  both.  If  geophysical  effects  caused  these  differences, 
the  mean  diurnal  amplitude  over  a  two-year  period  would 
have  to  be  greater  than  10  K  with  minima  occurring  near 
local  noon.  Recently  States  and  Gardner  [1998,  2000b]  ob¬ 
served  diurnal  amplitudes  of  this  order  at  these  altitudes 
during  most  seasons  but  with  temperature  maxima  occur¬ 
ring  at  local  noon.  Climatological  comparisons  of  data  from 
ground-based  lidars  and  daytime  Solar  Mesosphere  Explorer 
(SME)  observations  by  Clancy  et  al .  [1994]  also  show  that, 
on  average,  mid- latitude  temperatures  are  warmer  during 
daytime  at  87  km.  In  addition,  discontinuous  jumps  in  the 
monthly  average  local  time  behavior  of  the  three  data  sets 
are  observed  at  sunrise  and  sunset.  By  adding  bias  correc¬ 
tion  terms  to  the  inversion,  the  global  mean  of  each  data  set 
is  self-consistent ly  adjusted  so  that  the  local  time  behavior 
is  continuous  at  the  dawn/dusk  boundaries  between  data 
sets. 

Figure  2  shows  the  seasonal  evolution  of  the  diurnal  har¬ 
monic  estimated  with  and  without  a  bias  correction  term  in 
(1).  The  results  obtained  without  a  bias  correction  (Fig.  2a) 
indicate  a  diurnal  harmonic  as  large  as  15  K,  with  near  con¬ 
stant  phase  maxima  at  night.  The  results  obtained  with  the 
bias  correction  term  (Fig.  2b)  indicate  phase  maxima  oc- 
curing  mostly  during  the  day,  supporting  the  conclusions  of 
States  and  Gardner  [1998,  2000b]  and  Clancy  et  al.  [1994]. 
In  addition  the  winter  pre-dawn  temperature  enhancement 
in  Figure  2b  is  also  observed  by  States  and  Gardner  [2000b]. 
This  is  thought  to  be  the  effect  of  in-situ  chemical  heating 
from  the  H+O3  reaction  which  peaks  in  winter  at  night-time 
at  87  km  [Mlynczak  et  al .,  1991].  We  also  note  that  with  or 
without  a  bias  term,  our  estimates  of  the  seasonal  evolution 
of  the  semidiurnal  tidal  amplitudes  and  phases  are  generally 
consistent  with  other  estimates  in  the  literature  [  Williams 
et  al .,  1998;  States  and  Gardner ,  2000b]. 

The  derived  systematic  bias  correction  factors  defined  in 
(1)  are  -7.2  ±  3.4  K  (HRDI),  8.7  ±  4.7  K  (CSU)  and  10.5 
d=  4.7  K  (PMI).  In  reporting  these  values  we  have  assumed 
that  systematic  biases  are  equally  probable  in  each  measure¬ 
ment  technique.  The  sample  mean  and  standard  deviation 
of  the  total  data  set  is  195.3  K  and  19.7  K,  respectively. 
Based  on  the  magnitudes  of  instrument  biases,  the  system¬ 
atic  uncertainty  in  the  annual  mean  temperature  at  87  km 
over  this  period  is  approximately  =b  8.5  K.  Thus,  by  a  sta¬ 
tistical  consensus  of  currently  available  data  the  magnitude 
of  the  mid-latitude  annual  mean  temperature  at  mesopause 
altitudes  is  uncertain. 

Figure  3  shows  the  total  day- of-year /local- time  temper¬ 
ature  estimates  at  87  km  at  0°E  (Figure  3a)  and  180°  E 
(Figure  3b).  These  fields  have  been  reconstructed  using  the 
estimated  Fourier  coefficients  in  (1).  The  error-propagated 
uncertainties  of  the  individual  harmonic  components  are  on 
the  order  of  d=  3  K.  On  comparing  Figure  3a  and  Figure 
3b,  we  see  that  the  seasonal  variations  observed  at  two 
hypothetical  ground  stations  separated  by  180°  appear  to 
be  different.  Thus,  two  ground-based  observers  situated  at 
these  locations  would  infer  different  inter- annual,  annual  and 
semi-annual  harmonics  even  if  24  horn*  measurements  were 
available  to  remove  tidal  effects. 

The  reason  for  this  is  illustrated  in  Figure  4.  It  shows 
the  quasi-st  at  ionary  planetary  wave  field  reconstructed  us¬ 
ing  only  the  longitude/day-of-year  harmonic  terms  in  (1). 
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Figure  1.  The  location  of  the  satellite-  and  ground-based  obser¬ 
vations  in  the  space-time  (upper)  and  time-time  (lower)  domain. 
The  HRDI  observations  are  shown  in  blue,  CSU  lidar  in  red,  and 
PMI  in  green. 
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Figure  3.  Reconstruction  of  the  mid-latitude  87-km  day-of- 
year/local  time  temperature  surface  at  0°  E  (upper)  and  180°  E 
(lower)  for  1993  and  1994. 


During  late  1993  and  early  1994  significant  stationary  lon¬ 
gitudinal  variations  on  the  order  of  8  to  10  K  are  present. 

Smith  [1997]  observed  similar  planetary-scale  waves  at  these 
latitudes  and  heights  in  HRDI  wind  data  during  early  1994 
and  showed  that  they  are  strongly  correlated  to  stationary 
planetary  waves  observed  at  10  mb  (~  30  km).  Figure  4  also 
shows  that  stationary  planetary  waves  with  amplitudes  on 
the  order  of  5  K  occur  during  other  seasons.  The  existence 
of  these  stationary  planetary  waves  in  the  data  was  verified, 
independent  of  our  fitting  analysis,  by  simply  binning  the 
HRDI  observations  in  longitude  over  a  36-day  period.  We 
also  noted  that  longitudinal  temperature  differences  seen  in 
the  HRDI  data  between  the  CSU  and  PMI  sites  were  consis¬ 
tent  with  the  longitudinal  temperature  differences  observed 
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Figure  2.  The  seasonal  variation  of  the  diurnal  harmonics  in  Figure  4.  Reconstruction  of  the  mid-latitude  stationary  wave 
K  estimated  without  (upper  panel)  and  with  (lower  panel)  bias  at  87  km  (day-of-year  versus  longitude),  showing  significant 

correction  factors.  quasi- stationary  planetary  wave  variations. 
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in  the  ground  stations  after  removing  the  local  time  varia¬ 
tions  and  biases. 

The  annual  temperature  variation,  as  well  as  significant 
local  time  variations,  can  be  seen  in  the  temperature  fields 
in  Figure  3.  Due  to  the  superposition  of  tidal  modes,  and  in- 
situ  chemical  heating,  daily  local  time  fluctuations  on  the  or¬ 
der  of  zb  10  K  are  observed.  The  analysis  suggests  that  these 
variations,  longitudinal  asymmetries,  and  instrument  biases 
all  play  a  role  in  current  climatological  estimates  determined 
from  individual  ground-  or  space-based  instruments.  For  ex¬ 
ample,  after  self-consist ently  accounting  for  these  effects  in 
our  inversion,  the  amplitude  of  the  zonal  mean  annual  har¬ 
monic  at  87  km  is  estimated  to  be  14.1  it  2.1  K  with  a 
phase  maximum  occurring  at  356.2  ±6.2  days.  This  ampli¬ 
tude  is  less  than  17.0  ±  2.0  K  derived  by  She  and  vonZahn 
[1998]  using  nighttime  CSU  lidar  data  and  20.0  ±  0.2  K 
by  Niciejewski  and  Killeen  [1995]  using  PMI  data.  We  also 
note  that  this  amplitude  is  significantly  greater  than  the  8 
K  estimate  of  States  and  Gardner  [2000a]  made  during  the 
1996-1998  period  using  a  24  hour  lidar.  It  is  not  unreason¬ 
able  to  believe,  however,  that  the  amplitude  of  the  annual 
oscillation  at  87  km  could  vary  by  as  much  as  ±  5  K  on 
inter  annual  scales. 

4.  Conclusion 

To  produce  reliable  estimates  of  the  spatiotemporal  vari¬ 
ability  of  the  MLT  region  we  have  presented  a  methodology 
for  combining  satellite-  and  ground-based  data  sets.  It  was 
applied  to  determine  self-consist  ently  the  seasonal  evolution 
of  the  mean  state,  tides,  and  stationary  planetary  waves  at 
87  km  at  mid-latitudes  to  an  estimated  uncertainty  of  ±  3  K. 
Because  of  calibration  biases  between  the  available  satellite- 
and  ground-based  instruments,  we  were  unable  to  establish 
accurately  an  annual  mean  temperature  (193.5  K)  to  better 
than  ±  8.5  K.  However,  by  considering  the  relative  varia¬ 
tions  in  the  data  we  resolved  spatiotemporal  variations  in 
greater  detail  and  accuracy  than  would  have  been  possible 
using  data  from  a  single  instrument. 

The  derived  tidal  amplitudes  and  phases  support  the  re¬ 
sults  of  recent  observations  and  tidal  analyses.  Further,  our 
results  also  indicated  that  in-sit u  chemical  heating  may  be 
important.  We  also  observed  significant  stationary  plane¬ 
tary  waves  in  MLT  temperatures  (amplitudes  ~10  K).  The 
magnitude  of  these  previously  unreported  features  suggests 
that  they  play  an  important  role  in  the  dynamics  and  en¬ 
ergetics  of  the  region.  Furthermore,  these  waves  have  the 
potential  to  alias  into  climatological  estimates  of  seasonal 
fluctuations  made  from  single  ground  stations. 
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